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ON  THE  RESPONSE  OF  AIRPLANES  IN  A  3-DIMENSIONAL  GUST  FIELD 

by  John  C.  Houbolt 
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ABSTRACT 


The  response  of  an  airplane  in  a  3-dimensional  gust  field  is 
examined  with  the  intent  of  reducing  the  problem  to  the  simplest 
terms.  It  is  shown  that  response  evaluation  may  be  reduced  to 
separate  treatment  of  the  longitudinal  and  lateral  response  cases. 
The  longitudinal  case  reduces  to  considering  the  vertical  gust 
velocities  only,  and  generally  the  degrees  of  freedom  of  vertical 
motion,  pitch,  and  vertical  bending  modes;  vertical  motion  is  the 
prime  degree  of  freedom.  The  lateral  case  reduces  to  side  gust 
velocities  only,  with  degrees  of  freedom  of  side  motion,  yaw,  and 
side  bending;  yaw  is  the  prime  degree-of-freedom  in  this  case.  The 
longitudinal  gusts  are  found  to  be  unimportant.  The  analysis  shows 
that  the  results  from  the  separate  longitudinal  and  lateral  response 
evaluations  must  be  combined  to  obtain  the  fuselage  loads  and  loads 
on  T-type  tails. 


SYMBOLS 


a . 

av 

b 


c 


H(a>) 

L 

s 

S 


slope  of  the  lift  curve 

slope  of  the  lift  curve  associated  with  horizontal  tail 
slope  of  the  lift  curve  associated  with  vertical  tail 
span 
chord 

mean  aerodynamic  chord 
moment  arms 

frequency  response  functions 
lift,  also  turbulence  scale 
separation  distance 
wing  area 


SrjSv  areas  associated  with  horizontal  and  vertical  tails 
t  time 

u  longitudinal  component  of  turbulence 

U  airplane  velocity 

v  lateral  component  of  turbulence 

w  vertical  component  of  turbulence 

W  airplane  weight 

x,y,z  coordinate  axes  system;  also  displacements 

a  angle  of  attack 

T  dihedral  angle 

0  pitch  angle 

p  air  density 

ow  rms  value  for  vertical  gusts  w 

<p  roll  angle 

0  power  spectrum  of  vertical  gust  velocities 

w 

V  yaw  angle 

a)  circular  frequency 

Q  spatial  frequency,  ^ 


INTRODUCTION 


This  report  deals  with  the  response  of  an  airplane  in  a  3- 
aimensional  gust  field.  Various  aspects  of  the  problem  have  been 
considered  in  previous  studies,  references  1-17,  but  coverage  re¬ 
lative  to  the  1-dimensional  response  case  (vertical  gusts  random 
in  the  flight  direction  only)  is  small,  particularly  with  respect 
to  practical  applications.  The  studies  that  have  been  made  give 
consideration  to  such  items  as  the  derivation  of  cross-spectra  for 
isotropic  turbulence,  the  effect  of  considering  the  gusts  to  be 
random  in  the  spanwise  direction,  as  well  as  the  direction  of 


2 


and 


flight  (2-d  gust  field),  longitudinal  and  later al  response, 
horizontal  and  vertical  tail  loads. 

The  main  purpose  of  this  report  is  to  ^iv®  a  general  assess¬ 
ment  of  the  3-<3  turbulence  encounter  problem  with  the  aim  of  re¬ 
ducing  response  considerations  to  the  simplest  terms.  The  relative 
magnitudes  of  the  various  forcing  terms  due  to  each  of  the  three 
gust  components,  u  ,  v  ,  and  w  ,  are  developed.  The  terms  that 
are  felt  to  be  of  significance  are  singled  out.  Reduction  to  the 
subcases  of  longitudinal  and  lateral  response  is  made,  but  the 
essential  aspects  of  3*d  encounter  are  retained.  The  means  for 
treating  tail  and  fuselage  loads  due  to  the  combined  action  of 
vertical  and  horizontal  gusts  is  developed,  with  particular  atten¬ 
tion  being  givtin  to  T-taile,  wh?-€  a  superposition  of  vurtioal 
loadings  due  to  both  side  and  vertical  gusts  is  involved. 


THE  3-D  GUST  FIELD 


Figure  1  depicts  the  three -dimensions!  gust  field  that  is 
considered  in  this  analysis.  The  longitudinal  component  u  is 


considered  to  be  random  in  the  direction  of  flight  and  in  the  span- 
wise  direction.  This  component,  though  included  at  the  start,  is 
shown  to  be  of  negligible  importance  in  response  treatments.  The 
side  gusts  v  are  considered  to  be  random  in  the  flight  direction, 
but  due  to  the  smallness  of  the  airplane  in  depth,  are  considered 
uniform  in  the  vertical  direction.  The  vertical  component  w  is 
considered  random  in  both  the  flight  direction  and  in  the  spanwise 
direction.  Subsequently,  v  and ‘  w  components  are  shown  to  be 
the  prime  components  that  must  be  included  in  response  evaluation. 

The  gust  field  is  considered  to  be  isotropic.  Various  studies 
show  that  the  spectral  relations  due  to  von  Karman  apply  quite  well 
to  atmospheric  turbulence,  and  hence  they  are  adopted  as  the  basis 
for  the  mathematical  modelling  of  the  turbulence  field.  The  re¬ 
lations  are  (references  18  and  19): 


to*  w  and  v  , 

0w(fiL) 

for  u  , 

0^  (fZL) 


1  +  1  (1.339AL)2 

plll/b 

1  t  (1.339fiL) 


57S 


1  +  (I.339OL)2 


(1) 


(2) 


Isotropy  implies  that  the  severity  of  all  thr^e  components  is  equal 
thus,  as  expressed  by  the  rms  values 


u 


=  a  -  a. 


w 


(3) 


Isotropy  also  means  that  the  cross-correlation  and  cross-spectra 
between  the  three  components  vanish;  thus 


A  —  <b  —  (b  =0  (4) 

vuv  ^uw  yvw  K  ' 

These  relations  have  significant  consequence,  since  response 
evaluations  are  thereby  greatly  simplified. 


In  addition,  the  condition  of  isotropy  allows  for  the  con¬ 
venient  evaluation  of  the  cross-spectra  for  the  individual  gust 
components.  Two  of  the  more  important  cross-spectra  cases  are 


depicted  in  the  following  sketches 


In  case  (a)  we  are  interested  in  the  cross-spectra  between  the  w 
values  along  path  1  and  the  w  values  along  path  2,  as  would  be 
involved  in  the  consideration  of  nonuniform  spanwise  gusts.  Be¬ 
cause  of  isotropy,  we  can  derive  this  function  directly  from  the 
point  correlation  function  R(x)  .  We  replace  x  by  r  to  obtain 
the  cross-correlation  function 


The  cross-spectra  follows  as  the  Fourier  transform  of  this  function. 
Reference  16  applies  this  technique  to  derive  the  cross-spectral 
functions  associated  with  equations  (1)  and  (2). 

In  case  (b)  we  are  interested  in  the  cross-spectra  between  the 
w  values  at  point  1  and  the  w  values  at  point  2,  as  would  be  in¬ 
volved  in  considering  the  gust  loading  effects  on  the  wing  and 
horizontal  tail.  In  this  case,  it  is  common  to  consider  that  the 
w  values  at  point  2  are  simply  the  w  values  that  were  sensed  at 
point  3*  a  time  At  earlier,  where  At  is  the  time  required  for  a 
point  on  the  airplane  to  travel  the  distance  s  between  the  points, 

or  At  =  ~  .  (This  is  the  temporarily  frozen  gust  field  concept.) 

In  terms  of  the  point  correlation  function,  the  cross-correlation 
function  foi  this  case  is  simply 


5 


R12(t,s)  =  R(t  -  §) 


Prom  this  relation  the  cross-spectra  for  case  (b)  follows  as 


<*>12(a>,s)  =  e 


■jCUS 
"  * 


0(0)) 


where  0(to)  is  simply  the  point  spectrum,  equation  (1). 

Application  of  cross-spectra  in  multidimensional  response 
treatment  can  be  made  by  the  general  methods  outlined  in  reference 

1. 


GENERAL  RESPONSE  CONSIDERATIONS 


In  multidimensional  response  studies  wherein  the  airplane  is 
treated  as  a  rigid  body,  a  stability-axes  system  has  often  been 
used,  so  that  direct  use  of  stability  derivatives  could  be  made. 

We  shall,  however,  speak  mostly  in  terms  of  a  rectangular  coordinate 
system  x-y-z  moving  with  the  airplane  velocity  U  ,  but  otherwise 
inertially  fixed*  x  is  the  flight  direction,  y  the  spanwise 

direction,  and  z  is  vertical.  We  assume  that  if  information  on 

the  stability  derivatives  of  the  airplane  is  available,  it  can  be 
incorporated  by  a  suitable  axes  transformation. 

In  general,  the  airplane  is  considered  to  respond  in  three 
basic  ways 

a)  translation  in  the  x  ,  y  ,  and  z  directions 

b)  rotation  about  the  x  ,  y  ,  and  z  axes 

c)  flexible  body  deformations  . 

As  a  way  of  discussing  the  response  problem  in  general,  we  speak 
in  terms  of  the  basic  normal  mode  approach.  Thus,  the  following 
equation  applies 


M  a  +  p  a 
n  n  n  n 


A  a  =  f 
n  n  J 


pz  dS 
^  n 


This  equation  applies,  of  course,  to  all  six  rigid-body  modes,  as 
well  as  any  flexible  mode  that  is  included.  The  three  gust  com¬ 
ponents  enter  through  the  manner  in  which  they  give  rise  to  the 
pressure  term  p  that  appears  in  the  integral.  This  pressure 
term  also  includes  all  forces  that  develop  as  a  result  of  aircraft 
motion.  Through  use  of  equation  (5),  general  response  treatments 


may  be  made,  such  as  given  in  reference  1. 


RELATIVE  SIZE  OP  FORCING  FUNCTIONS 


In  this  section  we  consider  the  application  of  equation  (5)  to 
all  the  rigid-body  modes,  and  establish  how  each  of  the  three  gust 
component^  u  ,  v  ,  and  w  ,  influence  the  magnitude  of  the  forcing 
term  that  appears  on  the  right-hand  side. 

To  give  a  quick  overall  perspective  of  the  relative  magnitude 
of  the  various  forcing  terms,  we  give  at  the  onset  a  summary  of  the 
results  of  this  section.  Table  I.  This  table  indicates  the  order 
of  magnitude  of  the  forces  and  moments  that  are  applied  to  the  air¬ 
craft  due  to  each  of  the  three  gust  components;  W  is  aircraft 
weight,  b  the  span,  a  is  the  angle  of  attack  necessary  to  sus¬ 
tain  level  flight  of  the  airplane  at  the  speed  U  .  In  the  develop¬ 
ment  of  the  table  it  was  assumed  that  both  the  left  and  right  halves 
of  the  wing  experience  different  gust  velocities,  and  thus  the 
effects  of  spanwise  variations  in  turbulence  are  also  brought  out . 
Example  interpretations  are  as  follows.  The  magnitude  of  the  force 
in  the  z  direction  due  to  u  is  seen  to  be  essentially  the  weight 
times  a  characteristically  small  angle;  if  u^  and  u2  have  rms 
values  of  3  fps,  and  U  is  500  fps,  then  a  rms  force  in  the  z- 
direction  due  to  u  of  only  about  .01  W  is  indicated.  By  con¬ 
trast  consider  the  z-force  due  to  w  ;  we  see  the  force  is  given 
by  the  weight  times  the  ratio  of  comparable  angle  of  attacks.  Thus 
the  z-force  due  to  w  is  on  the  order  of  the  weight  of  the  air¬ 
plane.  Throughout  the  derivations,  the  weight  of  the  airplane  ex¬ 
pressed  in  terms  of  the  steady-state  lift  is  used  as  a  convenient 
reference  force  magnitude;  hence,  the  expression 

L  =  W  =  |  pU2Sa  (6) 


was  employed  where  a  is  slope  of  the  lift  curve,  S  is  wing  area 
and  a  is  the  steady-state  angle  of  attack  for  level  flight. 

Since  we  seek  only  to  establish  the  magnitude  of  the  forces  and 
moments  that  act  on  the  airplane,  without  being  precise  in  detail, 
we  use  an  analysis  akin  to  strip-theory,  and  assume  that  the  gusts 
act  in  quasi-steady  sense  only. 

Force  in  the  x-direction.-  It  is  usually  assumed  that  the  gust 
components  u  ,  v  ,  and  w  cause  negligible  perturbations  in  the 
x-direction,  and  this  assumption  is  also  adopted  here.  The  es¬ 
sential  reasoning  behind  this  quite  plausible  assumption  is  that 
the  airplane  in  the  x-direction  is  a  purposely  tailored  stream¬ 
lined  shape.  Of  the  three  components,  u  causes  the  largest 
variation  in  the  x-force,  or  drag,  but  the  force  is  small  as  the 


following  analysis  shows.  Consider  the  drag  of  the  airplane  ex¬ 
pressed  as 

D  =  D^U2 

With  the  use  of  differential  calculus  techniques,  incremental 
changes  in  drag  due  to  changes  in  U  are  given  by 


AD  =  2D1UAU 

If  we  combine  these  two  equations,  and  introduce  the  fact  that 
L  =  W  ,  we  find 

AD  D  2AU 

W  L  U 


If  AU  is  regarded  as  the  gust  component  u  ,  then  we  see  that  only 
small  drag  variations  in  terms  of  airplane  weight  are  involved;  for 
example,  a  rms  value  of  u  =  3  fps  ,  U  =  500  fps  ,  =  8  ,  indi¬ 
cates  D 

AD  =  .0015  W 

Since  this  drag  variation  is  small,  and  since  the  x-force  varia¬ 
tions  due  to  v  and  w  are  even  smaller,  we  arrive  at  the  con¬ 
clusion  that  we  can  ignore  the  x-force  due  to  all  three  components, 
as  the  table  indicates. 

Force  in  the  y-direction.-  No  fundamental  mechanism  exists  for 
producing  force  in  the  y-cirection  due  to  u  and  w  ;  therefore, 
a  blank  spot  for  both  of  tiese  terms  is  shown  in  the  table.  Tne 
v  component  can  produce  a  small  side  force  by  acting  on  the  fuse¬ 
lage  and  vertical  tail.  With  reference  to  the  following  sketch 


the  side  force  due  to  v  is  given  by 

av  p  y 

Fyv  -  -r  Pu  sv  § 


8 


where  a  and  Sy  are  associated  with  an  equivalent  vertical  tail 
which  produces  the  same  side  force  as  does  the  actual  fuselage  and 
tall  combination  of  the  airplane  when  in  an  angle  of  yaw  equal  to 


Z  .  With  the  use  of  equation  (6),  this  equation  becomes 

-  a  S  ,  - 

v  =  JL  JL  1 1  w 

yv  a  S  a  U 


(7) 


is  small;  for 
,  a  side 


This  relation  shows  th$t  the  side  force  due  to  v 

av  i  S  i  v 

example,  with  ~  2  9  y  ±y  j  and  even  with  ^  = 

force  equal  to  1/lb  the  weight  is  shown.  We  suggest  later  the 
omission  of  this  force  in  first-order  response  evaluations. 

Force  in  the  z-direction.-  The  force  in  the  z-direction  due  to 
u  isderived  as  follows.  The  lift  on  the  wing  is  given  by  (with 
the  use  of  strip  theory) 


L  = 


b/2 

fp  f 


cu12ady 


-b/2 


where  u^  is  the  instantaneous  flow  velocity  over  a  chord  due  to 
the  steady  forward  velocity  U  and  the  longitudinal  gusts  u  ; 
that  is 

u^  =  U  +  u 


If  this  expression  is  substituted  in  the  integral,  then  the  change 
in  L  from  the  steady-state  value  is  found  to  be  (ignoring  the 


u‘ 


term) 


u 


b/2 

2  p  J  c2Uuady 
-b/2 


With  equation  (6),  this  equation  may  be  written 

1 


'u 


=  „  fiLH 

y  =o u 


b/2. 


(8) 


where  c  is  the  mean  aerodynamic  chord.  To  bring  out  the  results 
a  little°more  clearly  we  apply  this  equation  to  the  idealized  wing 
shown  in  figure  2.  In  this  figure,  we  consider  u±  to  be  the  u 


p 

}j 

[j 

D 

0 


0 


L 


Figure  2.  Wing  system  used  to  derive  gust  forcing 
functions  -• 


velocity  that  acts  over  the  left  half  of  the  wing,  and  u2  the  u 
velocity  that  acts  over  the  right  half.  For  this  case,  equation 
(8)  yields 


(9) 


which  is  the  value  shown  in  the  table.  The  F,  force  due  to  u 
is  thus  seen  to  be  a  small  component  of  the  weight,  roughly  the 
weight  times  a  small  angle  (an  angle  of  attack  given  by  2  2i)  . 
We  will  eventually  neglect  this  force.  u 

The  gust  component  v  produces  no  force  in  the  z-direction 
and  so  no  entry  in  the  table  is  made. 


The  gust  component  w  is  the  strongest  Fz 
equation  for  lift  produced  by  w  is  given  as 

b/2 


producer.  The 


L  =  F. 


w 


-fp/ 

-b/2 


cU2  g  dy 


With  equation  (6)  this  equation  becomes 

1 

_n  r 

a  2  J 


w 


-1 


c 

c 


w 

U  Ab72 


(10) 


1C 


|V I'HilTTfr'  ifon  II  i  im(  raiWnrir-Tlirf 


In  application  to  our  simple  model  wing  of  figure  2,  this  equation 
yields 


-50*$ 


(ID 


which  is  the  result  given  in  Table  I.  Notice  the  magnitude  of  this 


term  by  comparison  with  the 


term  for 


the  value  is 


essentially  the  F*  term  for  u  divided  by  a  small  angle,  and  so 
a  force  an  order  of  magnitude  greater  is  found.  Since  the  force  is 
equal  to  the  weight  multiplied  by  the  ratio  of  angles  of  attack, 
which  are  roughly  equal  to  one  another,  the  force  is  essentially 


the  weight  of  the  airplane.  The 
cant  in  response  evaluations. 


due  to  w  is  thus  signifi- 


Rolling  moment . -  We  consider  the  rolling  moments  due  to  u 
and  w  first.  These  rolling  moment  forcing  terms  come  from  the 
same  lifting  forces  as  are  involved  in  deriving  equations  (8)  and 
(10) .  For  roll,  a  moment  arm  y  would  also  appear  in  the  inte¬ 
grals;  thus,  the  rolling  moments  may  be  shown  to  be 


“r  -W1  f  th) 


=  ai  / f:  §  tk  Kbfe) 


(13) 


These  equations  indicate  that  the  forcing  moments  from  the  left 
and  right  wings  try  to  cancel  one  another.  Indeed,  if  the  gusts 
u  and  w  were  uniform  in  the  spanwise  direction,  no  rolling 
power  would  be  created.  For  u  and  w  variable  in  the  spanwise 
direction,  some  rolling  moment  can  develop.  The  application  of 
equations  (12)  and  (13)  to  the  example  wing  of  figure  2  gives  the 
results 


-  Mb/Jl  _  ^ 


M  Wb  1/^1  w2'\ 

X  =  v  -  rr  ) 


.r  _ ...  I^il-I-  -IMI 


which  are  the  results  given  in  Table  I.  We  hold  back  further  dis¬ 
cussion  of  the  magnitudes  of  these  terms  until  the  rolling  moment 
due  to  v  is  established. 

The  forcing  rolling  moment  due  to  v  is  associated  with  a 
dihedral  effect.  The  local  angle  of  attack  due  to  a  side  velocity 
v  can  be  shown  to  be 


where  r  is  the  dihedral  angle.  Since  the  dihedral  angles  of  the 
right  and  left  wings  are  opposite  in  sign,  the  local  angle  of 
attack  and,  hence,  the  lift  will  be  opposite  also.  Additive  rolling 
moments  due  to  the  right  and  left  halves  are  therefore  produced. 

The  rolling  moment  due  to  v  is 

b/2 

=  f  P  2  f  cu2r  ~  y  dy 
o 

which  by  equation  (6)  may  be  expressed 


(16) 


If  v  is  assumed  constant  across  the  span,  the  application  of  this 
equation  to  the  example  wing  of  figure  2  gives 


m  Wb  r  v 

\  “TJu 


(17) 


as  is  shown  in  Table  I. 

.  We  now  compare  the  magnitudes  of  the  three  terms.  The  factor 
•j7  can  be  considered  as  a  common  moment  arm;  thus  we  can  compare 

the  force  producing  the  rolling  moments.  For  we  see  the 

force  is  essentially  the  weight  W  times  the  difference  of  small 
angles;  in  general,  therefore,  only  a  small  force  is  involved.  We 
can  thus  infer  that  this  rolling  moment  term  can  be  dropped,  and 
later  we  will  do  so.  The  term  is  seen  to  be  essentially  the 

Mp  term  divided  by  a  ;  thus,  even  though  the  difference  of  small 
Ku 


1 

% 

.3 
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angles  is  involved,  the  term  is  an  order  of  magnitude  greater  than 
the  term.  It  is  therefore  usually  large  enough  so  that  it 

should  be  considered.  For  MR  ,  the  factor  r  essentially  takes 

the  place  of  the  difference  of  small  angles;  because  the  term  in¬ 
volves  division  by  a  ,  however,  just  as  with  the  Mr  term,  it 

can  have  a  magnitude  as  great  the  Mr  term.  The  Mr  term  must 
therefore  generally  be  retained.  w  v 


Pitching  moment.-  Because  of  lift  variations  similar  to  those 
described  under  the  F„  term  due  to  u  ,  a  small  pitching  moment 
due  to  u  may  occur;  the  moment  is  quite  small,  however,  and  may 
be  neglected,  as  the  table  shows.  There  is  no  pitching  moment  due  to 
v  ,  and  so  there  is  no  entry  for  this  term  either.  The  main  pitch¬ 
ing  moment  is  caused  by  w  .  For  this  treatment  we  assume  that  the 
moment  is  primarily  due  to  a  vertical  tail  force  that  develops  on 
the  horizontal  tail  and  which  acts  at  some  moment  arm  e„  from  the 
c.g.  This  moment  is  given  by  n 


M  _  aH  Q  tt2  w  a 
Mp  2  P  SHU  U  H 


By  equation  (6)  this  value  becomes 


=  W 


SH  1  w 
a  T  aU 


'H 


(18) 


which  is  the  result  shown  in  the  table.  If  we  consider  that  a 
is  roughly  equal  to  |,  this  equation  shows  that  the  forcing 

pitching  moment  is  a  good  percentage  of  the  weight  W  times  a  sub¬ 
stantial  moment  arm  eH  .  Thus,  the  pitching  moment  due  to  w  may, 
in  general,  be  of  concern. 

Yawing  moment.-  The  yawing  moments  are  analogous  to  the  pitch¬ 
ing  moments.  The  yawing  moment  due  to  u  is  negligible,  and  w 
does  not  create  any.  The  yawing  moment  due  to  v  is  similar  in 
development  to  the  pitching  moment  due  to  w  .  Essentially,  we 
consider  the  side  force  Fy  due  to  v  that  develops  on  the  verti¬ 
cal  tall,  and  assume  it  to^act  through  a  moment  arm  about  the  c.g. 
of  ey  .  The  result  is 


a  S  , 

v  v  1  v 

a  S  a  U 


v 


(19) 


as  is  shown  in  the  table.  This  moment  is  found  to  be  of  signifi¬ 
cance  in  lateral  response  considerations. 
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Vertical  motion  modes.-  As  indicated  in  Table  I,  only  w  has 
practical  significance  in  exciting  the  normal  modes  of  the  airplane 
which  are  characterized  primarily  by  vertical  motion.  The  extent 
of  excitation  of  a  normal  mode  is  governed  by  a  generalized  force 
(right  side  of  equation  (5))  which  depends  on  the  modal  shape  being 
excited.  Since  this  exciting  force  varies  with  each  modal  shape, 
no  magnitude  evaluations  are  made  for  entry  into  the  table. 

Side  bending  modes.-  For  the  modes  which  are  characterized 
mainly  by  side  bending  motion  of  the  fuselage,  only  v  is  signifi¬ 
cant  as  a  forcing  function.  We  indicate  this  fact  by  a  check  mark 
in  Table  I,  analogous  to  the  check  mark  under  w  for  the  vertical 
motion  modes. 


SEPARATION  OF  LONGITUDINAL  AND  LATERAL  RESPONSE 


When  we  examine  Table  I,  we  note  that  an  orderly  grouping  of 
terms  appear,  and  recognize  that  the  longitudinal  and  lateral 
response  cases  may  be  separated,  thus  leading  to  significant 
simplifications  in  3-d  response  evaluations.  Before  discussing 
this  separation,  we  make  some  observations  of  interest  about  Table 
I.  We  note  that  of  all  the  degrees  of  freedom,  only  roll  is  in¬ 
fluenced  directly  by  all  three  gust  components.  In  spite  of  this 
fadt,  we  will  later  suppress  roll  as  a  degree  of  freedom,  since  it 
does  not  have  a  great  influence  on  load  response  evaluations.  The 
previous  section  indicated  that  the  terms  under  u  are  small  and 
so  we  will  now  drop  them.  Reference  6  verifies,  in  fact,  that  the 
u  forcing  term  for  roll  has  very  little  effect  on  the  roll  re¬ 
sponse.  With  the  elimination  of  the  x  degree  of  freedom,  and  the 
u  forcing  terms,  the  general  rigid-body  response  equations  would 
reduce  to  the  form 


y 

F 

yv 

z 

Fz 

zw 

c 

0 

— 

1  v  rw 

e 

Mp 

w 

tam  am 

my 

iv 

(20) 


In  this  symbolic  representation,  the  matrix  C  is  considered  to 
include  the  inertia  terms,  cross  coupling,  and  all  force  and  moment 
terms  that  are  associated  with  airplane  motion. 
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If  the  angular  motions  are  small,  then  the  longitudinal  and 
lateral  motions  become  essentially  uncoupled  and  we  may  split 
equation  (20)  into  two  separate  smaller  sets,  one  for  longitudinal 
response,  one  for  lateral;  we  discuss  these  sets  separately. 


Longitudinal  response.-  If  we  extract  the  longitudinal  response 
equations  from  equation  (£o)  and  introduce  the  flexible  response 
modes,  which  were  temporarily  suppressed,  we  find  the  following  set 
applying  to  longitudinal  response 


z 

F 

zw 

e 

Mp 

rw 

A 

• 

= 

• 

an 

On 

nw 

• 

• 

(21) 


We  notice  that  only  w  is  involved.  Generally,  the  equation 
applies  to  a  w-gust  field  which  is  both  random  in  the  direction 
of  flight  and  the  spanwise  direction.  The  matrix  A  contains  in¬ 
ertia  terms  and  aerodynamic  terms,  associated  with  airplane  motion, 
as  found  by  lifting  surfaces  methods,  for  example.  Various  simpli¬ 
fications  involve  considering  the  gusts  to  be  uniform  spanwise,  and 
of  suppressing  various  degrees  of  freedom. 


Of  all  the  degrees  of  freedom  in  equation  (21),  z  ,  the  rigid- 
body  vertical  translation,  is  the  most  significant.  Consideration 
of  z  alone,  as  is  often  done,  leads  to  first  order  results  for 
loads;  the  equation  for  this  case  appears  simply  as 


A-.Z  =  F  (22) 

w 

In  spite  of  the  fact  that  we  refer  to  this  equation  as  a  rigid-body 
equation,  it  should  be  noted  that  if  this  equation  is  used  for  air¬ 
planes  with  swept  wings,  care  should  be  taken  to  use  the  slope  of 
the  lift  curve  for  the  flexible  airplane  so  that  steady  state  wing 
bending  effects  are  taken  into  account. 

Lateral  re sponse . -  The  lateral  response  equations  that  are 
indicated  by  equation  (20),  together  with  an  appropriate  flexible 
mode  response,  are 
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We  notice  that  the  lateral  response  is,  in  general,  influenced  by 
both  v  and  w  .  In  equation  (23)  we  indicate  only  one  flexible 
mode,  since  it  is  felt  that  the  use  of  one  should  be  adequate  for 
most  all  practical  applications. 

If  motion  response  of  the  airplane  is  of  prime  consideration, 
then  the  degree  of  freedom  of  roll  (0)  should  be  included  in 
applications,  so  as  to  adequately  represent  the  dutch  roll  mode  of 
the  airplane.  The  roll  freedom  may  be  suppressed,  however,  if  tail 
and  fuselage  loads,  rather  than  overall  airplane  motion,  are  of 
chief  concern,  as  is  usually  the  case;  these  loads  are  not  affected 
greatly  by  roll  motion.  Yawing  motion  appears  to  be  the  biggest 
factor  associated  with  the  loads  induced.  The  lateral  response  case 
treated  in  reference  20  was  in  fact  examined  on  this  point.  The 
axes  systems  were  first  transformed  from  the  stability  system  to  the 
system  of  this  paper,  using  the  transformation  v  =  y  -  Utf  .  It  was 
found  that  the  results  obtained  for  the  case  of  y  =  0  ,  using  only 
V  ,  were  essentially  the  same  as  given  in  the  reference .  With  0 
suppressed,  the  lateral  response  equations  become 


We  note  that  only  v  is  now  involved. 

The  most  significant  degree  of  freedom  for  lateral  response 
appears  to  be  the  rigid-body  yaw  (^);  the  equation  for  this  case 
is  simply 

Bp^  =  My  (25) 

Use  of  this  equation  should  give  good  first  order  results  for 
lateral  loads.  We  wish  to  emphasize  the  contrasting  significance 


!  n 


5 


i  f 
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of  equations  (22)  and  (25).  For  longitudinal  response,  rigid-body 
vertical  motion  is  predominant;  this  result  may  be  attributed  to 
the  fact  that  large  aerodynamic  surfaces  are  exposed  to  the  gusts 
w  .  For  lateral  response,  rigid-body  yaw  is  the  most  significant; 
in  this  case  the  amount  of  aerodynamic  surfaces  exposed  to  v  is 
small,  but  the  surfaces  that  are  involved  can  create  powerful  yawing 
moments. 


FUSELAGE  LOADS 


We  have  seen  that  the  3-d  gust  encounter  problem  may  be 
separated  into  longitudinal  and  lateral  response  considerations. 

The  wing  loads,  for  design  purposes,  are  established  by  the  longi¬ 
tudinal  response  analysis,  using  equations  (21)  or  (22).  The  fuse¬ 
lage  loads  should,  however,  be  established  by  combining  the  longi¬ 
tudinal  and  lateral  responses.  Even  though  the  responses  are 
established  in  an  uncoupled  sense,  a  combining  of  loads  is  necessary 
to  arrive  at  design  loads.  The  following  analysis  illustrates  the 
procedure . 

Consider  a  fuselage  cross  section  and  the  moments  that  act  at 
this  cross  section,  as  shown  in  the  following  sketch 


M. 


The  moment  Mx  is  considered  to  be  the  result  of  encountering  the 
gusts  w  •  My  is  due  to  v  .  These  moments  may  be  given  as 


x  =/ w(T)hx(t  -  t)  dt 
-00 

t 

[y  v(t)hy(t  -  t)  dT 


(26) 

(27) 
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where  hx  is  the  x-moment  that  develops  due  to  a  unit-iropuD.se  w 
gust,  ana  hy.  is  the  y-moment  that  develops  due  to  a  unit-impulse 
v  gust.  (These  impulse  functions  are  referred  to  here  in  implicit 
form  only,  since  they  are  really  not  derived.) 


The  stress  at  point  A  is  given  by 


s 


+ 


M 
-2L. 
I 


If  we  substitute  equations  (26)  and  (27)  into  this  expression,  and 
take  the  Fourier  transform,  we  obtain 


Fs  =  t  Vw  +  f-y  Vv 


where  Hx  and  Hy  are  the  frequency  response  functions  for  Mx 
and  My  due,  respectively,  to  unit  sinusoidal  w  and  v  gusts, 
and  and  Fv  are  the  Fourier  transforms  of  the  w  and  v 

fusts:  the  Hx  and  Hv  functions  are  the  important  response 
unctions  that  required  valuation.  The  spectrum  for  s  "follows 
from  this  equation  as 

2  2  2 
Hxl  ^w  +  lHyl  ^v 
y 


<t>Qn 


2 
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In  the  derivation  of  this  equation,  use  was  made  of  the  fact  that 
the  cross-spectrum  between  v  and  w  vanishes,  equation  (4).  Since 
0V  =  0W  ,  equation  (28)  may  be  written 

05(<d)  =  | h|  2  0W  (29) 

where 


The  rms  value  a  of  stress,  and  N0  value,  for  use  in  design, 
follow  directly  from  equation  (29)  in  the  customary  fashion. 

It  is  of  interest  to  note  how  the  spectral  results  for  stress 
differ  from  the  result  that  is  obtained  when  Mx  and  My  are 
static-type  moments.  Consider  the  cross  section  to  be  circular 
and  let 


x  =  r  cos  0 


y  =  r  sin  0 


Then,  if  the  moments  axe  static,  the  stress  is  given  by 


Mr  M  r 

s  =  sin  0  +  cos  0 


The  angle  leading  to  maximum  stress  is  given  by 


,  Mr 
ds  _  x 


cos  0  - 


M  r 

sin  0 

y 


We  solve  this  equation  for  0  ,  and  substitute  the  result  in  the 
equation  for  s  to  find  the  maximum  stress.  The  point  is,  when 
static  moments  are  involved,  there  is  a  position  0  which  leads  to 
maximum  stress,  as  might  be  expected. 

By  contrast,  the  spectral  result  for  stress  for  this  circular 
cross-section  case,  as  found  through  equation  (28),  is  of  the  form 

2  2  2  2  2 
ag  =  a1  sin  0  +  a2  cos  6 

To  find  the  location  of  maximum  a  we  differentiate  and  obtain 

s 


do  2  2 

■jg —  =  2(a^  “  °2  )  s*n  e  cos  6  =  0 

From  this  expression  we  see  that  maximum  cg  is  either  at  6=0 
or  0  =  ^  .  An  exception  occurs  when  ai  =  °2  >  ^is  case>  a£ 

is  invariant  with  respect  to  6  . 


TAIL  LOADS 


For  3“d  gust  response,  tail  loads  merit  a  special  treatment. 

For  conventional  t’ils,  no  unusual  problem  arises.  The  longitudinal 
response  involving  i  is  used  to  calculate  the  loads  on  the  hori¬ 
zontal  tail.  In  similar  fashion,  the  lateral  response  involving 
v  is  used  for  the  vertical  tall  loads.  For  T-tails,  load  determi¬ 
nation  is  more  involved  and  requires  the  combining  of  loads  due  to 


the  longitudinal  and  lateral  response.  The  following  analysis  de¬ 
picts  the  situation. 

Consider  a  T-tail  as  shown  in  the  following  sketch 


1*1  ^2 


and  consider  the  tail  to  have  a  pitch  attitude  0  and  a  yaw  atti¬ 
tude  V'  .  The  use  of  aerodynamic  interference  theory  indicates 
that  the  loads  shown  are  given  "by  equations  of  the  type 

=  a-^0  -  b-j^ 

L2  =  ax0  +  b-jtf  (30) 

L3  =  c -jtf 

Thus,  a  yawing  motion  creates  loads  on  the  horizontal  tail,  as  does 
motion  in  pitch.  It  should  be  mentioned  that  0  may  arise  from 
at  least  three  sources,  an  actual  pitch  position  0  ,  an  effective 

angle  of  attack  ^  due  to  vertical  motion,  and  a  gust  angle  of 

attack  ^  ;  similar  comments  apply  to  y  . 

By  equation  (30),  we  see  that  the  loads  on  the  horizontal  tail 
are  due  to  both  the  longitudinal  and  lateral  responses.  A  combining 
of  loads  must  therefore  be  considered ;  this  combining  process  is 
similar  to  that  given  in  the  previous  section  on  fuselage  loads. 

We  illustrate  it  here  in  connection  with  determining  the  bending 
moment  at  point  A  shown  in  the  sketch.  The  moment  is  given  by 

t  t 

MA  =  I  w(T)hw(t  "  T)dT  +J  v(T)hv(t  '  T)  dT 
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Reduction  of  this  equation  to  spectral  form  gives  the  result 

*MA  =  |Hj2  *w  +  |HV|2  *v 

where  Hw  and  Hy  are  the  frequency  response  functions  for  moment 
at  A  due  to  unit  sinusoidal  w  and  v  gusts,  respectively.  In 
the  derivation  of  this  equation  use  was  again  made  of  the  fact  that 
0Vw  =  0  ,  equation  (4).  Since  0V  =  ^  ,  the  equation  becomes 

%“Ohj2  +  ihvI2X  (3d 


Design  rms  values  for  bending  moment,  and  associated  N  ,  follow 
from  this  equation.  ° 

The  establishment  of  the  moment  at  point  B  does  not  involve 
w  .  Care  must  be  taken,  however,  to  include  the  load  Lo  ,  and  the 
loads  Lj  arid  L2  ,  when  deriving  this  moment;  in  this  instance, 
lateral  response  alone  is  involved. 


CONCLUDING  REMARKS 


This  report  considers  the  3- dimensional  gust  response  problem 
of  an  airplane.  The  aim  was  mainly  to  suggest  ideas  on  how  the 
problem  could  be  reduced  to  simplest  terms.  Table  II  is  presented 
by  way  of  summarizing  the  basic  notions  of  the  paper.  Some  of  the 
chief  points  are  as  follows: 

a)  The  x  degree  of  freedom  can  be  ignored. 

b)  The  longitudinal  gusts  u  may  be  ignored. 

c)  Longitudinal  and  lateral  motion  may  be  separated  and 
treated  in  an  uncoupled  way. 

d)  Only  the  vertical  component  w  enters  into  the  longi¬ 
tudinal  response  problem. 

e)  The  complete  lateral  response  which  includes  the  dutch 
roll  mode  involves  both  v  and  w  . 

f)  For  loads  determination  in  the  lateral  response  case, 
roll  may  be  suppressed,  and  hence  treatment  is  in  terms 
of  v  only . 
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g)  The  longitudinal  response  determines  the  wing  loads  for 
use  in  design;  for  the  very  large  aircraft,  consideration 
should  be  given  to  the  spanwise  variations  in  gusts, 
especially  if  the  turbulence  scale  is  found  to  be  less 
than  1000  ft. 

h)  Tall  and  fuselage  loads  are  found  by  combining  the  results 
that  are  obtained  from  separate  longitudinal  and  lateral 
response  evaluations. 

i)  First-order  load  results  can  be  obtained  by  considering 
rigid-body  vertical  motion  only  in  connection  with  w  , 
and  rigid-body  yawing  motion  only  in  connection  with  v  . 


TABLE  I 

MAGNITUDE  OP  FORCING  TERMS  DUE  TO  3-D  GUSTS 

u 

V 

w 

x  -  Force 

0 

0 

0 

y  -  Force 

0 

av  Sv  1  V 

W  JTB”  w  tr 

0 

z  -  Force 

■ft**) 

0 

W  l/Wl  W2^ 

?  3crr  +  tv 

Rolling  Moment 

«R 

* 

Wb^l  U2>| 

t \ir  ‘  tv 

Wb  r  v 
T"S  IT 

Wb  ]/wl  w2'\ 

t  ’  tv 

Pitching  Moment 

i  Mp 

0 

0 

w  aH  SH  1  w 

"  r  t  ^  u  6h 

Yawing  Moment 

0 

u  av  Sv  1  V 

"  a"T"3tr 

e  v  0 

Vertical  Modes 

0 

0 

V 

Side  Bending  Modes 

0 

V 

0 
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